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Abstract 

O . 

QT 1 We study perfect fluid Stephani quantum cosmological model. In the present work 

5^ ' the Schutz's variational formalism which recovers the notion of time is applied. This 

gives rise to Wheeler-DeWitt equation for the scale factor. We use the eigenfunctions in 
order to construct wave packets for each case. We study the time-dependent behavior 
of the expectation value of the scale factor, using many-worlds and deBroglie-Bohm 
• interpretations of quantum mechanics. 

(N ' 
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On 
O . 

O ! 1 Introduction 

> 

^ ■ In recent years observations show that the expansion of the Universe is accelerating in the 
present epoch [1] contrary to Friedmann-Robertson- Walker (FRW) cosmological models, with 
non-relativistic matter and radiation. Some different physical scenarios using exotic form of 
matter have been suggested to resolve this problem [2, 3, 4, 5, 6]. In fact the presence of exotic 
matter is not necessary to drive an accelerated expansion. Instead we can relax the assumption 
of the homogeneity of space, leaving the isotropy with respect to one point. The most general 
class of non-static, perfect fluid solutions of Einsteins equations that are conformally flat is 
known as the "Stephani Universe" [7, 8]. This model can be embedded in a five-dimensional flat 
pseudo-Euclidean space, which is not expansion-free and has non- vanishing density [9, 7, 10]. 
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In general, it has no symmetry at all, although its three dimensional spatial sections are 
homogeneous and isotropic [11]. The spherically symmetric Stephani Universes and some of 
their subcases have been examined in numerous papers [8]. So it may be important to study 
the quantum behavior of this models. 

The notion of time can be recovered in some cases of quantum cosmology, for example when 
gravity is coupled to a perfect fluid [12, 13, 14]. This kind of systems are often studied as follows 
[16, 13, 14]. First one uses the Schutzs formalism for the description of the perfect fluid [17, 18], 
second one selects the dynamical variable of perfect fluid as the reference time. Finally, one uses 
canonical quantization to obtain the Wheeler-Dewitt (WD) equation in minisuperspace, which 
is a Schrodinger-like equation [12]. After solving the equation, one can construct wave packets 
from the resulting modes. The wave packets can be used to compute the time-dependent 
behavior of the scale factor. If the selected time variable results in a close correspondence 
between the expectation value of the scale factor and the classical prediction (prediction of 
General Relativity) for long enough time, the selected time variable can be considered as 
acceptable. This approach has been extensively employed in the literature, indicating in general 
the suppression of the initial singularity [15, 16, 12, 19, 20, 13, 21, 14]. 

We can also study this situation from ontological interpretation of quantum mechanics 
[13, 22, 23, 24, 25, 26]. In this approach, the problem of time is solved in any situation (not 
only in the presence of the matter field). In fact, the ontological interpretation predicts that 
the system follows a real trajectory given by 

P q = S, q , (1) 

where p q is the momentum conjugated to the variable q. Here, S is the phase of the wave 
function (\l> = Re lS ), where R and S are real functions. The equation of motion (1) is governed 
not only by a classical potential V but also by the following quantum potential Vq = ■ 

In the present paper, we use the formalism of quantum cosmology in order to quantize 
Stephani cosmological model in Schutz's formalism [17, 18] and find WD equation in minisu- 
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perspace. In the Schutz's variational formalism the wave function depends on the scale factor 
a and on the canonical variable associated to the fluid, which plays the role of time T. Here, 
we describe matter as a perfect fluid matter. Although, this is semiclassical from the start, 
but it has the advantage of defining a variable, connected with the matter degrees of freedom 
which can be identified with time and leads to a well-defined Hilbert space. 

2 The Model 

The action for gravity plus perfect fluid in Schutz's formalism is written as 

S=[ d 4 x^jR + 2 [ d i xVhh ab K ab + ( d 4 xy/^p, (2) 

JM JdM JM 

where K ab is the extrinsic curvature and h a b is the induced metric over the three-dimensional 
spatial hypersurface, which is the boundary dM of the four dimensional manifold M in units 
where IQnG = 1 [27]. The last term of (2) represents the matter contribution to the total 
action, where p is the pressure. In Schutz's formalism [17, 18] the fluid's four-velocity is 
expressed in terms of five potentials e, (, £, 9 and r: 

Uv = -(e, u + (^ + 6r u ), (3) 

where p is the specific enthalpy, the variable r is the specific entropy, while the potentials ( 
and £ are connected with rotation and are absent in models of the FRW type. The variables e 
and 9 have no clear physical meaning. The four-velocity satisfies the following normalization 
condition 

u v u v = -1. (4) 

The metric in spherically symmetric Stephani Universe [9, 28, 7, 11, 8, 29] has the following 
form, 



ds 2 = - 



F(t) 



a{t) d [V(r,ty 



V(r,t) dt \ a(r,t) 



(5) 
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where the functions V(r,t) and F(t) are defined as 

1 



V(r,t) = l + - A k(t)r\ (6) 

F(t) = . a{t) (7) 
IC 2 (t)a 2 (t) - k(t) 



Using the line element (5) and the Einstein's equation, one can easily show the functions C(t), 
k(t) and a(t) are not all independent but are related to each other with the expressions 

3C 2 (t) 

M = ^ ( § ) 
p(r,t) = J-[2C(t)C(t) KMZ^l -3C 2 (t)], (9) 

where an overdot denotes a derivative with respect to t. Note that in the spherically symmetric 
Stephani models with the given coordinate system, the energy density pit) is uniform, contrary 
to the pressure p(r,t), which is not uniform and depends on the distance from the symmetry 
center at r = 0. This is the reason for the absence of the barotropic equation of state (p = p{p)) 
in such models. However, if we assume some relations between p(t) and p(r, t), this could allow 
us to eliminate one of the unknown functions such as C(t). Therefore we are left with two 
unknown functions k(t) and a(t). The first one k(t), plays the role of a spatial curvature, and 
the second one a(t), is the Stephani version of the FRW scale factor. 

Now, we consider an observer which is placed at the symmetry center of the spherically 
symmetric Stephani Universe and our physical assumptions will concern in the vicinity of 
r pa 0. First, we assume that locally, matter in the Universe satisfies the barotropic equation 
of state 

p(r Pa 0,t) = ap(t). (10) 

By substituting the Stephani metric (Eq. (5)) in the action (Eq. (2)) and Choosing a curvature 
function k(t) in the form [30] 

k{t) = (3a?{t), (11) 
and after some thermodynamical considerations [12], the final reduced effective action near 



r ?y 0, takes the form 



S = [ dt\ - 6^ + + iV" 1 ^ ^T7TT(^ + ^) 1/a+1 ex pf " ^ 



(12) 



(a + 1) 1 

The reduced action may be further simplified by canonical methods [12] to the super-Hamiltonian 



n 



„2 n a+l t 



(13) 



24a a J 

where p a = —12aa/N and p e = —pou°Na 3 , po being the rest mass density of the fluid. Using 
the canonical transformations 



T = p T e- T p; (a+1) , 
e = e — [a + 1) — , 

Pe 



PT = P : +1 e T , 



Pe = Pe: 



(14) 



which are the generalization of the ones used in Ref. [12], the super-Hamiltonian takes the 
form 
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Ti = -Ph - 6/3a 1+7 + , (15) 



24a a 3a 

where the momentum p? is the only remaining canonical variable associated with matter which 
appears linearly in the super-Hamiltonian. 

The classical dynamics is governed by the Hamilton equations, derived from Eq. (15) and 
Poisson brackets, namely 

d = {a,NH}-- 



Np a 
T2^ 



p a = { Pa ,NH} 



-^l + 6iV(l + 7)^ + ^, 



t = {T, NH} = Na~ 3a 
p T = { Pt ,NH} = 0. 



(16) 



We also have the constraint equation H — 0. Choosing the gauge N = a 3a , we have the 
following equations for the system 



T 



(17) 



a = (3a-h--l(l + 1 )Pa^- 1 -ja^- 2 p T , (18) 
I a I 4 

Note that the classical equations for 7 = +1 case in Ref. [31], are corresponding with choosing 
the gauge N — 1. In this gauge the constraint equation TC = reduces to 

-6aa 2 - 6/3a 7+1 + a~ 3 > T = 0, (20) 



or 



Imposing the standard quantization conditions on the canonical variables (p a — > — 2^-, — > 
—i-ffp) an d demanding that the super-Hamiltonian operator (15) annihilate the wave function, 
we are led to the following WD equation in minisuperspace (fi = 1) 

^— - 144/3a 2+7 tf - i24a 1 " 3a ^ = 0. (22) 
da 1 at 

According to the equation (17) T = t can be associated with the time coordinate [32, 33]. 
Equation (22) takes the form of a Schrodinger equation id^/dt = H^f. As discussed in 
[20, 33], in order for the Hamiltonian operator H to be self-adjoint the inner product of any 
two wave functions $ and \& must take the form 

roc 

= / a 1_3a $**da. (23) 
J 

On the other hand, the wave functions should satisfy the following boundary conditions [33, 34] 

d$(a,t) 



^(0,t) = or 



da 

The WD equation (22) can be solved by separation of variables as follows, 



= 0. (24) 

a=0 



q(a,t) =e iEt iP(a), (25) 
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where the scale factor dependent part of the wave function (tp(a)) satisfies 



ip" + 144/3a 2+ > = 24£a 1 " 3 >, 



(26) 



and the prime denotes derivative with respect to a. 

An interesting feature of the Stephani model is that the spatial curvature is time-dependent. 
The recent observational data show that our Universe is spatially flat. Moreover, negative pow- 
ers in equation (11) lead to the spatially flat Universe in the present epoch. In the next section 
we solve a general class of exactly solvable models which 7 = — (1 + 3a) and find the corre- 
sponding eigenfunctions. Then we construct the wave packets by appropriate superimposing of 
these eigenfunctions and compute the expectation value of the scale factor versus time. After 
solving the classical equations exactly, we compare the classical and quantum solutions and 
show that these solutions are asymptotically the same. Moreover, we study the case using 
de-Broglie Bohm interpretation of quantum mechanics and find the corresponding bohmian 
trajectories which are different from the classical case for small times due to the effect of the 
quantum potential. 

3 Quantum cosmological models with 7 = — (1 + 3a) 

In this Section, we consider the particular relation between 7 and a as 7 = —(1 + 3a). In this 
case the WD equation (26) takes the form 



and hence equation (25) has the following general time-dependent solutions under the form of 
Bessel functions 



ij" + 24(£ - 6p)a 



l-3a 



V^ = o, 



(27) 



3(l-a) 




96(£-6/3) 3(1-0) 




(28) 
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Wave packets can be constructed by superimposing these solutions to obtain physically allowed 
wave functions. The general structure of these wave packets are 

V(a,t)= A{E)ty E {a,t)dE. (29) 

J6/3 



We choose C2 = to satisfy the first boundary condition of Eq. (24). Defining r = ^ 9 ^f_ , 
simple analytical expressions for the wave packet can be found by choosing A(E) to be a 
quasi-gaussian function 

q,(a,t) = f ^^(^(H'+Wju^)^ (30) 
J o 

where v = 3 ^ l _ a ^ and k is an arbitrary positive constant. The above integral is known [35], 
and the wave packet takes the form 



3(l-a) 

-+6i/3t 



*(a,t)=a- AB 4 , 3Q , (31) 

(_2£)5(t=st 

where B = k — — a) 2 t. Now, we can verify what these quantum models predict for the 
behavior of the scale factor of the Universe. By adopting the many- worlds interpretation [36], 
and with regards to the inner product relation (23), the expectation value of the scale factor 

j^a l - 3a ^{a,tya^{a,t)da 



is easily computed, leading to 



F ( 3c *- 5 \ 18(l-a) 4 /2 , o (f 2 n _ 
I 3(a-l) ' r t "T *rx> \ 3(1 



(32) 



ft 



-p / 3a-4 \ 

1 U<*-i)J 

Now, we can calculate the dispersion of the wave packets 



3(l-a) 



(33) 



(Aa)? = (a 2 ). - (a) 2 , (34) 

using (31,32), we have 



This shows the dispersion of the wave packets through the time with the minimum at t — 0. 
This is similar to the free particle case, where the the wave packets disperse more rapidly for 
more localized initial states. 

The important feature of this model is the avoidance of the singularity. Equation (33) 
shows that the expectation value of the scale factor never vanishes for all time. On the other 
hand, at the quantum level, since the probability density of finding the scale factor at a = 
(with regards to the inner product relation (23) and the behavior of Bessel functions for small 
values of the argument) is zero in all times (lim a ^ a 1_3a |\l/(a, t)\ 2 = 0), we have a indication 
that these models may not have singularities at the quantum level. 

In classical case, by eliminating the pr variable in the equations of motions (18,19), the 
resulting equation in case 7 = — (1 + 3a) takes the following simple form 

a = (3a-l)^, (36) 



which has the exact solution as 
where 



a(t) = a t^-^, (37) 
r (^f)) rl8(l-a) 4 ' 

-p ( 3a-4 \ 
1 \3(a-l)J 



(32) 2 K 



(38) 



Figure 1 shows the behavior of the classical scale factor (37) and quantum mechanical expec- 
tation value of the scale factor (37) versus time for various cases. 

It is known that the results obtained by using the many-worlds interpretation agree with 
those that can be obtained using the ontological interpretation of quantum mechanics [22, 23, 
24, 25, 26]. In Bohmian interpretation the wave function is written as 

^ = Re iS (39) 

where R and S are real functions. Inserting this expression in the WD equation (22) (for 
7 = —(1 + 3a)), we have 




Figure 1: Classical scale factor and the expectation value of the scale factor versus time for 
a = 2/3 (up, left), a — 1/3 (up, right), a = (middle, left), a = —1/3 (middle, right), 
a = — 2/3 (down, left), and a = — 1 (down, right). 



10 



dR 



1 OR OS 



-R 



d 2 S 



0. 



(41) 



dt 12a 1 - 3 " da da 24a 1 ~ 3a da 2 
1 " R is the quantum potential which modifies the Hamilton- Jacobi equation. 



Here, Q — 2 Aa 1 -'^ r da 2 

In fact, the deviation from the classical trajectory happens whenever the quantum potential 
is more important than the classical potential. The real functions R(a, t) and S(a, t) can be 
obtained from the wave function (31) as 



R = 



S = 



v 



6(l-a) 



a exp 



na 



3(l-a) 



3 (1 - a) 2 a^- a H (4 -3a) 

+ irrz r arctan 



^ 2 +(|) (l-a)H 2 
3 (1 -aft 



128 



k 2 + 



1 - a) 4 t J 



3(1 - a) 



32 k 



+6/3*. 



(42) 
(43) 



In the Bohmian interpretation, the behavior of the scale factor is governed by the following 
equation 



OS 
da 



(44) 



On the other hand, from the definition of p a (16) and for N = a 3a we have p a = —12a 1 3a a. 



Therefore, Bohmian trajectory becomes 



512- = 3(1 -af 



K 



2 +(lJ (I"")** 5 



(45) 



which the integration yields 



a(t) = ao 



K 2 + 



32/ 



1 - a)H 2 



3(l-a) 



(46) 



where a is the constant of the integration. This completely coincides with the computation 
of the expectation value of the scale factor. Now we can find the quantum potential 
k 1 — a 



Q(a,t) 



32 



K 2 + 



[l-a)H< 



3k(1 - a)a 3 ^ - (4 - 3a) ^+(^) (I ~ a)H 2 }.(47) 



Using the relation between the scale factor and time (46) the quantum potential can be written 
in terms of the scale factor as 



Q(a) = k 



1 — a 3 (i_ Q ) (4 — 3a) — 3/t(l — a)a 



3(l-a) 




32 



a 3(l-a) 



(48) 
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It is obvious that the quantum effects are negligible for large values of the scale factor and are 
important for small values of the scale factor. Therefore, asymptotically we have the classical 
behavior. 

In the next section for completeness we consider briefly some interesting and exactly solvable 
cases which have bound state solutions. 

4 Bound state solutions 

In this Section, we study four different cases of 7 and a. We find the exact discrete energy 
spectrum and corresponding eigenf unctions. 

For 7 = — 1 and a = 1/3 (radiation), the WD equation (26) reduces to 



(49) 



which can be rewritten as 




(50) 



by taking x = a — 



E 



we have 



6J3 




(51) 



which is the Airy's differential equation [37]. This equation has two solutions as Ai ( 144/5) 1//3 x 
and Bi [( 144/5) 1//3 x . First one is exponentially decreasing function of x and the second one 
grows exponentially and is physically unacceptable. Therefore, the solution is 




(52) 



We choose the first boundary condition (24), which leads to 




(53) 
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. 4 
0.2 

to 

^ 
-0.2 
-0.4 

0.5 1 1.5 2 2.5 3 3.5 

a 

Figure 2: Plot of the wave function {ip(a)) for (3 = 1 and n = 8, showing the oscillatory 
behavior for the small values of the scale factor and exponential damping for the large values 
of the scale factor. 

Airy's function Ai(x) has infinitely many negative zeros z n = —a n , where a n > 0, therefore, 
the energy levels quantize and take the values 




The time-dependent eigenfunctions take the form 

tt n (a,t) = e iE "*Ai 

It is important to note that Airy's function Ai(x) has an oscillatory behavior for x < (a < |r^) 

E 

whiles for x > (a > -^k) decreases monotonically and is an exponentially damped function 
for large x (Fig. 2). Therefore, the solutions (55) show a classical behavior for small a and 
a quantum behavior for large a. This is contrary to usually expected results in the previous 
Section. In fact detecting quantum gravitational effects in large Universes is noticeable which 
has been also observed in FRW and Kaluza-Klein models [32, 38]. 

In 7 = — 1 and a = —1/3 (cosmic strings) case, the WD equation (26) reduces to 

-if," + lUpaip = 2AEa 2 i>. (56) 

For E < the above equation can be written as 

-^" + 24|£ 





3/f 
~E 




^ = 0, 



(57) 
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38 

by taking x = a — jt- we have 



d 2 21 6B 2 

; 4>(x) + 2A\E\x 2 4>{x) = f^f-^(x). (58) 



This equation is identical to the time-independent Schrodinger equation for a harmonic oscil- 
lator with unit mass and energy A: 

d 2 ijj{x) 22 



<ix 2 



+ wV^(x) = 2A^(a;), (59) 



216/3 2 

where 2A = and w 2 = 24|_E|. Therefore, the allowed values of A are uj(n + 1/2) and the 



possible values of E are 



E "-{7m^Y)) • n=0,1,2 (60) 

therefore, the stationary solutions are 

* B (a, = e*^Vn (a - ^ , (61) 

y? n (a;) = tf n ((24|£|)^) ^v 7 ^* 2 , (62) 

where if n are Hermite polynomials. The wave functions (61) are similar to the stationary 
quantum wormholes as defined in [39]. However, neither of the boundary conditions (24) can 
be satisfied by these wave functions. 

For 7 = — 2 and a = (dust regime), the WD equation (26) reduces to 

-ip" + lU(3ip = 24Eaip, (63) 

which can be written as 

V>" - 2AE (f -aU = 0, (64) 



by taking i=^-awe have 



< - 2AEx^ = 0, (65) 
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which is again the Airy's differential equation [37]. Therefore, the physically acceptable solution 
is 



ip(a) = Ai 



(24E7) 



1/3 W 



We choose the first boundary condition (24), which leads to 



Ai 



0. 



(66) 



(67) 



Airy's function Ai(x) has infinitely many negative zeros z n , therefore, the energy levels quantize 
and take the values of 



(66 



which exist only for negative values of f3. The time- dependent eigenfunctions take the form 



#„(a,f) = e l * nt Ai 



,1/3 f'M 

Ik 



— a 



(69) 



(24E n ) J 

These solutions (69) show a quantum behavior for small a and a classical behavior for large a. 
In 7 = — 2 and a = — 1/3 (cosmic strings) case, the WD equation (26) reduces to 

+ 144/3^ = 24Ea 2 ilj. (70) 

For E < the above equation can be written as 

-i>" + 24|£|aV> = -144/3^. (71) 

This equation is identical to the time-independent Schrodinger equation for a harmonic oscil- 
lator with unit mass and energy A, where 2A = —144/3 and uo 2 = 24|£'|. Therefore, the allowed 

values of A are u(n + 1/2) and the possible values of E are 

216/5 2 



E„ 



(n + 



l\2 



n = 0, 1,2, 



(72) 



for (3 < 0. Thus the stationary solutions are 

M/„(a,t)=e^Vn(«), (73) 

<p n (a) = H n ((U\E n \)*a) e~^ a \ (74) 

The solutions for odd n satisfy the first boundary condition (24) and the appropriate wave 
packets can be constructed by superposing these stationary solutions. 
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5 Conclusion 



In this work we have investigated perfect fluid Stephani quantum cosmological models. The use 
of Schutz's formalism allowed us to obtain a Schrodinger-like WD equation in which the only 
remaining matter degree of freedom plays the role of time. We have obtained eigenfunctions 
and therefore constructed the acceptable wave packets by appropriate linear combination of 
these eigenfunctions. The time evolution of the expectation value of the scale factor has been 
determined using the many-worlds and Bohmian interpretations of quantum cosmology. We 
have shown that contrary to the classical case, the expectation values of the scale factor avoid 
singularity in the quantum case. At the end, we solved some interesting bound state cases 
and found their discrete energy eigenvalues. We have also shown that in some bound state 
cases, we may observe the quantum effects in the large scales which correspond to the quantum 
behavior at the late time cosmology. 
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